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With advances in fabrication technologies it is now possible to create precisely controlled geome- 
tries and pinning landscapes for vortex matter in type-II superconductors. Here we use numerical 
simulations to examine vortex states and dynamics in periodic funnel geometries where a drive is 
applied in the easy flow direction. We show that this system exhibits a number of different com- 
mensurability effects when the vortex configurations match to both the periodicity of the array and 
the geometry of the funnels. The vortex configurations in this system are generally different from 
those observed for single isolated triangular superconducting samples due to the coupling of vortices 
in adjacent funnels. At certain matching fields, peaks in the critical current are absent due to the 
particular vortex configurations that occur at these fields. We find that the overall depinning force 
increases with increasing vortex density as a result of the enhanced vortex- vortex interactions caused 
by a crowding effect at the funnel tips. When a system becomes less mobile as a result of increased 
particle interactions, it is said to exhibit a jamming behavior. Under an applied drive we observe a 
series of elastic and plastic vortex flow phases which produce pronounced features such as jumps or 
dips in the transport curves. In all of the flow phases, only one vortex can pass through the funnel 
tip at a time due to the vortex-vortex repulsion forces. As a consequence of this constraint, we 
observe the remarkable result that the sum of the vortex velocities at a fixed drive remains nearly 
constant with increasing magnetic field B rather than increasing linearly. This result is similar to 
the behavior of sand in an hourglass. We also show how noise fluctuations can be used to distinguish 
the different flow phases. Our results should be readily generalizable to other systems of particles 
flowing in periodic funnel geometries, such as colloids or Wigner crystals. 

PACS numbers: 74.25.Uv,74.25.Wx 



I. INTRODUCTION 

Advanced nanostructuring techniques permit the cre- 
ation of specific superconducting structures for control- 
ling vortex motion in type-II superconductors. Experi- 
ments and simulations for square and triangular artifi- 
cial pinning arrays demonstrate that pronounced com- 
mensurability effects occur when the number of vortices 
is an integer multiple of the number of pinning sites, 
and that these effects can be observed as peaks in the 
critical current. 1-6 Above the first matching field, mul- 
tiple vortices may occupy individual pinning sites^ or 
only a portion of the vortices may be captured by the 
pinning sites with the remaining vortices located in in- 
terstitial sites where they can still be pinned by the re- 
pulsive interactions with vortices at the pinning sites 
Vortex matter in these periodic pinning arrays can also 
exhibit a remarkably rich variety of distinct dynami- 
cal phases originally predicted in simulations^^ These 
include one-dimensional flow of interstitial vortices be- 
tween immobile vortices in the pinning sites, disordered 
or turbulent phases where the number of moving vor- 
tices fluctuates strongly, and laminar states where the 
vortex flow can organize along the rows of pinning sites. 
The transitions between these different states appear as 
specific features in the voltage-current curves, includ- 
ing negative differential conductivity where the num- 
ber of moving vortices or the average vortex velocity 
decreases with increasing drive. In very recent experi- 
ments, these dynamical phases were observed in both low 
temperature^ and high temperature 12 superconductors 



with periodic pinning arrays. It is also possible to create 
periodic pinning arrays that contain intrinsic asymmetry, 
such as with asymmetric thickness modulation ; 13 ' 14 fun- 
nel geometries ^r— composite pinning sites^ or arrays 
of triangular traps ii2r— This asymmetry can produce a 
diode effect when the depinning force is higher in one 
direction, and can give a ratchet effect in which a net 
dc vortex flow occurs upon application of an ac drivel 
Reversals of the ratchet flow from the easy asymmetry 
direction to the hard asymmetry direction can occur as 
a function of magnetic field and other par ameter a 14 ' 18 ~ ^£ 
due to various collective interactions of the vortices. One 
of the earliest proposals for a vortex ratchet involved a 
periodic asymmetric channel or funnel^ and such ge- 
ometries have now been experimentally fabricated . 16 ' 17 ' 24 
The first experiments on this system verified the exis- 
tence of ratchet effects, but also detected boundary ef- 
fects caused by the disordered injection of vortices at the 
edge of the sampled The edge effect can be overcome by 
using sample geometries where the vortices flow in annu- 
lar asymmetric channels so that vortices need not enter 
or exit the sampled Such a technique may also provide 
the resolution required to investigate individual vortex 
motion through the channels^ 

In this work we study the dynamics of vortices in a 
periodic asymmetric channel geometry or funnel array 
and show that a variety of new types of vortex dynam- 
ics and behaviors can arise, including a jamming effect 
and flow patterns that are organized such that only one 
vortex at a time passes through the funnel tip. Such 
a geometry could be realized by etching the edge of a 
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single superconducting strip into a periodic funnel shape 
so that the vortices would flow in a true single chan- 
nel. Funnel geometries created with periodic channels 
in two-dimensional superconducting samples should also 
exhibit many of the same properties we observe; however, 
in these systems vortices located between the asymmet- 
ric channels could be important and under high drives 
could depin. The presence of vortices in the regions be- 
tween the channels can be deduced from the onset of 
hysteresis in the critical current curves I C (H) at fields 
higher than the fields at which many of the ratchet ef- 
fects are observed* 1 ^ Although our study focuses on in- 
teracting vortices in type-II superconductors, we expect 
that the same dynamics will occur for other systems 
of repulsively interacting particles in funnel geometries. 
These include charged colloidal assemblies ; 26 ! 27 magnetic 
colloids, 28 charged metallic dots,— and classical electron 
crystals. 30 Additionally, in recent experiments of ion flow 
through a single funnel, effects such as negative differ- 
ential conductivity were observed^ Our results suggest 
that the same type of clogging effect we observe could be 
occurring in these artificial ion channels. 



II. JAMMING AND CLOGGING IN VORTEX 
MATTER 

An important feature that differentiates funnel geome- 
tries from the other vortex ratchet geometries is that in 
the funnel the vortices are forced to move through a nar- 
row bottleneck. At this constricted point, the repulsive 
vortex-vortex interactions are very important and favor 
the organization of the vortex flow into a pattern that 
permits only a single vortex to pass through the bottle- 
neck at a time. The bottleneck has many features in 
common with granular hopper geometries where grains 
flow through a funnel. In the granular case, it is known 
that decreasing the width of the hopper aperture can 
cause the flow of grains to be impeded or jammedj 32 i 33 
Systems that become immobile due to particle-particle 
interactions are often referred to as jammed. The physics 
of jamming has attracted growing interest as a way to un- 
derstand many types of loose particle assemblies such as 
grains, colloids, and emulsions in situations where these 
systems exhibit a sudden onset of resistance to shear, 
with possible connections to the glass transitioni 34 i 35 In 
a granular hopper where grains flow through a thin fun- 
nel, the jamming effect occurs when some grains block 
the motion of other grains. The interaction between the 
grains is short ranged and has a sharp cutoff. Vortices 
also experience a repulsive interaction; however, it is sig- 
nificantly longer in range and smoother than the inter- 
action among grains. Understanding how systems with 
intermediate range or longer range interactions can de- 
velop a jammed state is an open question. 

Generally, in a vortex system, increasing the effective 
vortex-vortex interaction strength reduces the effective- 
ness of the pinning and causes the depinning force to 



decrease. If vortex matter can exhibit a jamming effect, 
the opposite behavior would occur and the system would 
become more immobile with increasing vortex- vortex in- 
teraction strength. We note that in the peak effect phe- 
nomenon, the effective pinning force increases with in- 
creasing vortex density or temperature. Many explana- 
tions of this effect involve the reduction of the effective 
vortex- vortex interaction force with increasing density or 
temperature due to changes in the penetration depth or 
softening of the vortex lattice near H C 2 or in order 
to match the normal expectation of increased pinning 
force caused by decreased vortex interaction strength. 
Superficially, however, the increased pinning force gener- 
ated by increased vortex density resembles a jamming ef- 
fect. In simulations of vortex systems with more vortices 
than pinning sites; 3 - 7 - gradually increasing the strength 
of the vortex-vortex interactions initially increased the 
depinning force since the vortices at the pinning sites 
blocked the free vortices from moving, similar to a jam- 
ming effect. On the other hand, when there are more 
pinning sites than vortices, the depinning force monoton- 
ically decreases as the vortex-vortex interaction strength 
increases^ 6 - This shows that if vortices are to exhibit jam- 
ming behavior, it must arise from the collective interac- 
tions of the vortices rather than from vortex-pin interac- 
tions alone. In our system we fix the bare vortex- vortex 
interaction strength and study the appearance of jam- 
ming and clogging effects due to density-induced changes 
in the effective importance of the vortex- vortex interac- 
tions. 



III. SIMULATION TECHNIQUES 

Following previous techniques for simulating vortices 
in periodic pinning geometries, we employ Langevin dy- 
namics in a two-dimensional system^ ' 10 i 14 ' 15 i 18 i 22 i 38 We 
consider a sample containing a single channel composed 
of N c = 16 funnels. Each funnel has a small aperture 
size of a = 1.8A, a wide aperture size of b — 7 AX, and a 
length of L c — 9A, where lengths are measured in units of 
the London penetration depth A. The funnels are aligned 
in the x direction and the sample has periodic boundary 
conditions along the x direction only. A total of N v vor- 
tices are placed only inside the funnel channel and the 
region outside the channel is empty. The motion of the 
vortices is calculated by integrating the following over- 
damped equation of motion: 

dR Nl ' 

'7^ 1 = -E VC/ ™(^)~ F wall+FD+Ff. (1) 

Here the damping constant is rj = 0od/27r£ 2 pjv in a crys- 
tal of thickness d, where 4>q = h/2e is the elementary flux 
quantum, £ is the superconducting coherence length, and 
Pn is the normal state resistivity. The vortex- vortex force 
is repulsive with a potential U vv (Rij) = AqK^R^ /A), 
where Aq = </>g/27r/ioA 3 , Kq is a modified Bessel function, 
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R-i(j) is the position of vortex and Rij — |Rj— Rj|. If 
a vortex approaches one of the channel walls sufficiently 
closely, it experiences a wall force F wa n. The channel 
walls are constructed out of Nj, = 4N C repulsive elon- 
gated potential barriers which are inverted versions of 
the potential wells employed in Ref. 39 . Each barrier has 
a central rectangular region which repels vortices in the 
direction transverse to the long direction of the rectan- 
gle, along with two half-parabolic cap regions which repel 
vortices from the ends of the barrier. We have 

N b 

K,n = (f P /r P ) £ Rf k e(r p - R% )6(4 - «R± 

k 

+ Rt k Q(r p - Ri)e(l k - i4)H*k- (2) 

Herei?± = |R,-R£±/ fc pf |, R^ = KRi-R^-p^l, R p k 
is the position of the center point of barrier fc, r p = 0.4A 
is the barrier radius or half width, f p = 15/o is the barrier 
strength, l k is half the length of the central rectangular 
region of barrier fc, and p| (p^ ) is a unit vector parallel 
(perpendicular) to the axis of barrier fc. The individual 
barriers are connected together to form a pair of sawtooth 
shapes as illustrated in the inset of Fig. [TJ The barrier 
lengths are 21 ^ = 2.8A for the vertical walls and 21^ — 
18A/-\/3 for the slanted funnel walls. The central portion 
of the channel is featureless and the vortices experience 
confining forces only from the barrier walls and not from 
the channel itself. The barriers are sufficiently strong 
that vortices can never cross them under the conditions 
considered in this work. 

The initial vortex positions are obtained by plac- 
ing the vortices evenly throughout the funnel and per- 
forming simulated annealing. Temperature is modeled 
as Langevin kicks F T with the following properties: 
(F T (t)) = and (F t T (t)FT(t')) = 2ryfc s T% ( 5(t - f), 
where ks is the Boltzmann constant. After the vortex 
positions are initialized, we apply an external drive F^ = 
FjjH representing the Lorentz force from an applied cur- 
rent in the positive x or easy-flow direction and measure 
the average vortex velocity (14) = N^ 1 Y^f v v « ' x - The 
drive is slowly increased in small increments with a fixed 
waiting time between each increment. The waiting time 
is taken sufficiently long that the system always reaches 
a steady state velocity at each drive before we make our 
measurements. The depinning force F c is defined as the 
drive at which (V x ) > 0.001. In this work we examine 
only the dynamics in the easy flow direction and note 
that in general F c is higher for driving in the hard or 
negative x direction so that a diode effect is possible, 
in agreement with experiments ! 16 ! 24 For driving in the 
hard direction, clogging type dynamics occur which will 
be explored elsewhere.— 

In our system the penetration depth A is less than the 
funnel size, so many of our results should carry over to 
typical colloidal systems with optical trap arrays. In the 
vortex ratchet experiments of Refj^, the size of the tri- 
angular traps is less than A so we are working in a differ- 
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FIG. 1: The depinning force F c vs B/B^, for vortices in a pe- 
riodic funnel array for driving in the +x, easy-flow direction. 
A series of peaks appear at matching fields of B/B^ = 4, 5, 
7, 8, 11, 12 and 14. Weaker peaks occur at B/B,/, = 3, 10, 
and 13, while peaks are absent at B/B^, = 2, 6, and 9. Inset: 
A portion of the sample showing the funnel array geometry. 



ent regime from the experiment, although we expect that 
many of the same effects can appear in both systems. Ex- 
periments with thicker films or larger funnel length scales 
should also be possible which would be much closer to the 
regime we consider. 



IV. COMMENSURATION EFFECTS IN THE 
EASY FLOW DIRECTION 

In Fig. [T] we plot the the depinning force F c for driv- 
ing in the easy +x direction vs B / B^, where B^ denotes 
the field at which there is one vortex per funnel. In our 
geometry, F c = at B/B$ = 1.0 since each vortex can 
move unimpeded through the bottleneck and then flow 
freely along the center axis of the funnel. Fig. [1] shows 
peaks in F c at B/B^ = 4.0, 5.0, 7.0, 8.0, 11.0, 12.0, 
and 14.0. Smaller peaks in F c appear at BjB^ = 3.0, 
10.0, and 13.0, while peaks are absent at BjB^ = 2.0, 
6.0, and 9.0. The varying sizes and shapes of the corn- 
mensuration peaks are consistent with the existence of 
distinct arrangements of singly quantized vortices within 
the funnel plaquettes at different matching fields. 

Commensurability effects in two-dimensional square or 
triangular periodic pinning arrays can take two differ- 
ent forms depending on whether multiple vortices are 
trapped at each pinning site or whether interstitial vor- 
tices are present. If multiple vortex pinning occurs, then 
the vortex configuration at each matching field B/B^ — 
n, with integer n, is the same as the configuration at 
B/Bj, — 1.0, but with rt-quantized vortices trapped at 
each pinning site. The result is a peak in the critical 
current at every matching field. 41 Similarly, the vortex 
configurations at fractional fields such as B/B^, = 1/2 
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are repeated at all fields B/B^ — n + Such multi- 

quanta commensuration effects are very similar to the 
commensuration effects observed in superconducting wire 
networks 42 If the pinning sites can capture a maximum 
of one vortex, then for fields above B/B^ = 1, some 
vortices will be located in the interstitial regions and the 
vortex lattice structures can be different at each matching 
field. 3,4 At matching fields where the interstitial vortex 
structure is disordered, there is no peak in the critical 
current^,. The magnitude and shape of the commensura- 
tion peaks show striking variations when different types 
of vortex crystals form at different matching fields. For 
example, in square pinning arrays, a square vortex lattice 
forms at B/B^ = 2.0, a less stable dimer lattice with a 
smaller critical current peak appears at B/B,f, = 3.0, a 
very stable triangular lattice with a strong critical cur- 
rent peak is present at B/B^ — 4.0, and the partially 
disordered vortex structures at B/B^ = 6.0 and 7.0 pro- 
duce no peaks in the critical current £ 

For the asymmetric funnel array, the quasi-one- 
dimensional nature of the system might be expected to 
produce identical commensuration effects at each match- 
ing field; however, it is possible for the vortices within 
each funnel to distort in both the x and y directions in or- 
der to try to form triangular ordering on a local scale, and 
thus the response of the system differs at different match- 
ing fields. Simulations and experiments on single meso- 
scopic triangular superconducting samples have shown 
that the vortices can form triangular or partially ordered 
configurations at magic fillings such as B/B^ = 3.0, 6.0, 
and 10.0.— At these matching fields, we find weak or 
missing commensuration peaks in the funnel geometry, 
as indicated in Fig. [1] The high symmetry of these well 
ordered states results in poor pinning of the vortices in 
the easy-flow direction since it allows some vortices to si- 
multaneously align along the x axis of the channel while 
closely approaching the narrow aperture of the funnel. 
The alignment of the vortices produces an additional re- 
direction force on the vortex closest to the funnel tip, 
permitting it to flow out of the funnel at a relatively 
low driving force. In contrast, at commensurate fillings 
where the vortices adopt partially disordered configura- 
tions or have degenerate ground states, the critical cur- 
rent is high. This is likely due to the lack of a well-defined 
easy shear direction for the disordered vortices. 

Figure [1] shows a vanishing critical current at low fields 
which results when the vortices are able to sit along the 
center of the funnel channel and can flow along the chan- 
nel unimpeded by the funnel geometry. The experiments 
of Refjii revealed a finite critical depinning force at low 
fields. Since the experimental channels are less than A 
wide, this finite critical force could be the result of a 
bowing effect in the center of the funnel channel that 
creates some effective pinning even for very low fields. 
The edge barrier may also be playing a role at low fields, 
and it is likely that some intrinsic random pinning exists 
throughout the entire sample which would create a finite 
F c at all fields. If the intrinsic pinning is weak, it should 




FIG. 2: Vortex positions (dots) and funnel geometry (lines) 
in a small portion of the sample at different fields with no 
applied drive, (a) At B/B^ — 2.0, an aligned dimer state 
forms, (b) At BjB$ — 3.0, a triangular state forms with 
the triangle orientation alternating in every other funnel, (c) 
B/B,p = 4.0. (d) B/B^ = 5.0. (e) At B/B^ = 6.0, a trian- 
gular structure forms that fits well with the triangular funnel 
geometry and that matches the vortex configuration found 
in an individual mesoscopic triangular superconductor. This 
field also corresponds to a missing commensuration peak in 
Fig- ED (f) B/B,,, = 7.0. 



be possible to observe the commensuration effects shown 
in Fig. [T] The intrinsic pinning effects can be strongly 
suppressed near T c , and in a later section we show that 
the commensuration effects in Fig. Q] are robust at finite 
temperature. We also note that for colloidal particles in 
an optical funnel trap array, such intrinsic pinning effects 
would not be present. 

At incommensurate fields n < B / B^ < n + 1, due to 
the discreteness of the vortices the funnels are occupied 
by a mixture of n and n+1 vortices. The depinning force 
is reduced at these fields as a result of the asymmetric 
repulsive force experienced by a vortex at the boundary 
between two funnels of different occupancy. Fig. [1] indi- 
cates that peaks or enhancements of the depinning force 
can arise at non-matching fields such as at B/B^ = 3.5 
and B/B^ — 5.75; however, unlike the matching peaks, 
these non-matching peaks are not robust against thermal 
fluctuations. 

In Fig.[2|a-f) we illustrate the vortex configurations at 
zero applied drive for B/B^ — 2.0, 3.0, 4.0, 5.0, 6.0, and 
7.0. For B/B^ = 2.0, Fig. Ufa) shows that each funnel 
captures two vortices which form a dimer state. All of 
the dimers are aligned in the same direction with one 
vortex in each dimer located at the upper left corner of 
the funnel and the other vortex located along the lower 
wall. We note that this is not the same configuration 
that would arise for a single isolated triangular super- 
conductor, where the two vortices would maximize their 
spacing by sitting with one vortex at the tip of the tri- 
angle and the other vortex in the center of the opposing 
triangle wall. Interactions between vortices in adjacent 
funnels in our system would make such an arrangement 
energetically unfavorable since the vortex at the tip of 
the triangle would be too close to the leftmost vortex in 
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the adjacent funnel. The tilted dimer configuration in 
Fig. HJa) minimizes the vortex-vortex interactions both 
within a single funnel and in adjacent funnels. 

The vortex dimer state in Fig.[^a) is two- fold degener- 
ate; in the other possible orientation, a vortex is located 
at the lower left corner of each funnel. In an infinitely 
long system at finite temperature or in the presence of 
quenched disorder, it is possible that domain wall excita- 
tions could form where the dimer orientation flips from 
one ground state to the other. In this case, it may be 
possible to map the dimer state to a one-dimensional 
Ising model which is known to have a long-range or- 
dered ground state only at T — 0.0. The mapping of 
effective dimer and trimer states of particles in periodic 
substrates to Ising and other spin models has been pro- 
posed previously for colloids on two-dimensional periodic 
substrates^ 4 - and vortices in honeycomb pinning arrays. 45 
For BjBfy = 3.0, Fig. [^b) shows that the three vortices 
in each funnel form a triangle with one vortex located in 
the corner of the funnel. The orientation of the triangle 
alternates in every other plaquette from having the upper 
funnel corner occupied by a vortex to having the lower 
funnel corner occupied by a vortex. This ground state 
has similarities to a one-dimensional antiferromagnetic 
ordering. 

For B/B^ = 4.0, shown in Fig. Htc), both corners of 
each funnel are occupied by vortices and the ground state 
is non-degenerate. This is the first filling at which a 
pronounced peak in F c emerges, as seen in Fig. [TJ The 
ground state at B/B^, = 5.0, illustrated in Fig. [5{d), is 
very similar to the configuration at B/B^ = 4.0 with 
the addition of one vortex in the center of the chan- 
nel near the tip of the funnel, while at B/B^ = 6.0, 
shown in Fig. HJe), there are two vortices in the center 
of the channel. At B/B^ — 7.0, Fig. [^f ) indicates that 
the configuration changes from states with four vortices 
along the walls and the remaining vortices in the cen- 
ter of the channel to a state with six vortices lining the 
walls and only one vortex in the center of the channel. 
Figure [1] indicates that there is a pronounced peak in 
F c at B/B<t> = 7.0 but not at B/B^ = 6.0. The vortex 
configurations at B/B^ = 2.0, 3.0, 4.0, 5.0, and 7.0 dif- 
fer from the configurations found in isolated mesoscopic 
triangular superconductors j 4 ^ however, the configuration 
at B/B^ = 6.0 is almost the same as that in an isolated 
triangle since the vortices can form an almost perfect tri- 
angular ordering within the funnel at this field. Since the 
accommodation of the triangular ordering to the bound- 
aries is so energetically favorable, it overcomes the energy 
cost of placing two vortices close together near the funnel 
aperture with one vortex shifted slightly in the positive 
x-direction. Since this vortex experiences an extra force 
from the other vortices, the effective pinning potential at 
the tip of the funnel for this vortex is depressed, lowering 
F c . At B/B^ = 7.0 this condition no longer holds since 
the extra vortex that had been near the funnel tip shifts 
to a new location along the wall, where it no longer exerts 
an extra force on the vortex at the center of the channel. 
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FIG. 3: Vortex positions (dots) and funnel geometry (lines) 
in a small portion of the sample at different fields with no 
applied drive, (a) B/B^ = 8.0. (b) B/B^ = 9.0, where there 
is a missing peak in F c . (c) B/B^ = 10.0. (d) B/B^ = 11.0. 
(e) B/Bcf, = 12.0. (f) B/B^ = 13.0. 



In Fig. [3ja), the ordered vortex configuration at 
B/Bfj, = 8.0 has four vortices on the upper wall of each 
funnel and three vortices on the bottom wall, along with 
one vortex in the center of the channel near the wide end 
of the funnel. This configuration has a two-fold degen- 
erate ground state since either the upper or the lower 
funnel wall could be occupied with the four vortices. At 
B/B<p = 9.0, Fig.^b) shows that there are now four vor- 
tices lining both the top and bottom funnel walls. The 
vortices on one wall are slightly more compressed than 
the vortices on the other wall so that the vortices near 
the funnel tip are not aligned in the y-direction. The 
wall with the stronger compression alternates from top 
to bottom in adjacent funnels. Figure [1] indicates that 
there is no peak in F c at B/B^ = 9.0. The relatively 
low depinning threshold at this field is a result of the 
close proximity of a vortex near each funnel tip to the 
vortex in the center of the channel combined with the 
asymmetry of the vortex compression along the walls; 
the extra force experienced by the channel vortex from 
the less compressed wall of the funnel causes it to de- 
pin more readily. At BjB^ = 10, shown in Fig. [3Jc), 
two vortices occupy the center of the channel and the 
remaining eight vortices line the upper and lower funnel 
walls in a symmetric configuration. In Fig. [3f d) we plot 
the configurations at B/B^ = 11.0, where four vortices 
line each funnel wall and a triangular vortex structure 
forms near the wide end of the funnel. At B/B^ = 12, 
Fig. ISJe) shows that an asymmetric configuration of five 
vortices on one funnel wall and four vortices on the other 
forms along with a skewed triangle of vortices in the open 
portion of the channel. At B/B^ = 13, the symmetric 
configuration illustrated in Fig.[3jf) forms with five vor- 
tices along each funnel wall and a triangle of vortices in 
the center of the channel. 
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FIG. 4: The average distance to the closest neighbor (d m in) 
at Fd = F c vs B/B^. The line indicates a fit to 1/B. 



A. Jamming and Bottleneck Effects Due to 
Vortex- Vortex Interactions 

Figure [T] shows that the depinning force tends to in- 
crease overall with increasing B. This is the result of a 
crowding effect that makes it more difficult for the vor- 
tex structure to distort in order to permit individual vor- 
tices to pass through the bottleneck. We can show that 
the crowding is geometrically induced by measuring the 
average closest neighbor distance (dmin) for the vortex 
configuration at the depinning force Fd = F c . To deter- 
mine (d m in), we solve the all- nearest- neighbors problem 
with a simple algorithm. The distance from each vortex 
i to its closest neighbor at R^ n is d min = |Rj — R^J. 

Then, (d min ) = N^ 1 ^f" d l min . Tne value 01 (dmin) at 
the depinning threshold has a convincing 1/B form, as 
shown in Fig. [4] As expected in a system where the 
pinning originates from vortex- vortex interaction forces, 
at a given field B a vortex can depin when the driving 
force pushes it closer to its neighboring vortex than the 
average spacing between vortices. As B increases, the av- 
erage vortex spacing should drop as 1/B, consistent with 
the behavior of (d m i n ). The critical force curve shown in 
Fig. [T] is very far from being a smooth function of B, un- 
like (d m in), and this simply indicates that the particular 
geometric arrangements of the vortices at different fields 
may make it easier or more difficult for two vortices to 
approach each other closely enough to depin. 

For the most commonly studied types of pinning, such 
as random pinning or arrays of individual pinning sites, 
when the number of vortices exceeds the number of pins, 
the depinning force tends to decrease with increasing B 
due to the relative increase in the strength of the vortex- 
vortex interactions compared to the pinning energy. In 
this the vortex lattice becomes stiffer, some vor- 

tices are forced to shift out of the pinning sites and oc- 
cupy interstitial regions. In other words, a stiff vortex 



lattice cannot adjust to the pinning site configuration as 
well as a soft vortex lattice can. One proposed mecha- 
nism for the peak effect observed in superconductors with 
random pinning is a softening of the vortex lattice due 
either to thermal fluctuations or to changes in A which 
increase the effectiveness of the pinning. In the funnel ge- 
ometry we consider here, depinning does not require the 
vortices to overcome the pinning strength of individual 
pinning sites. Instead, the vortices depin once they are 
able to overcome the vortex-vortex interactions blocking 
the passage of individual vortices through the funnel tips. 

A system that becomes less mobile due to in- 
creased interactions among the particles is said to be 
jammedj 34 ! 35 ' 46 ! 47 Studies on systems with short range in- 
teractions where a single probe particle is pushed through 
a collection of other particles in the absence of pinning 
have shown that the probe particle motion at a finite and 
constant driving force F^ becomes slower for increasing 
particle density, and at a critical density the probe par- 
ticle becomes stuck or jammedi 26 ' 35 i 48 In this jammed 
state, a critical driving force F c must be applied to unjam 
the probe particle, and this critical force monotonically 
increases with increasing density^ 8 - The critical force is 
analogous to the depinning force needed to move the vor- 
tices through the funnels in our system, which also shifts 
to higher values as the vortex density increases. There 
are many similarities between Fig. [T]and the behavior of 
jamming in colloidal systems. At very low densities, the 
critical depinning force drops to zero since the vortices 
are so far apart that they interact only extremely weakly. 
Previous studies of probe particles in systems with longer 
range particle-particle interactions found that a finite de- 
pinning force, similar to a jamming effect, exists even at 
low densities and increases monotonically with particle 
density. A key ingredient for jamming in our system is 
that a portion of the vortices must remain immobile. The 
vortices lining the walls of the funnel are held in place 
due to the strength of the repulsive vortex- vortex inter- 
actions; when one of these vortices moves toward the tip 
of the funnel, it experiences a barrier due to the compres- 
sion of the vortices at the funnel tip, providing a finite 
depinning force. Vortices away from the funnel walls are 
also pinned by means of this compressive repulsion. If 
the funnel walls did not converge, but instead remained 
a fixed distance apart, the depinning force would be ab- 
sent. In experiments and simulations on straight chan- 
nels where there are vortices both inside and outside the 
channels, the vortices outside of the channel are strongly 
pinned and create a periodic potential modulation for the 
vortices within the channel which allow the channel vor- 
tices to be pinned^ For the funnel geometry we consider, 
additional immobile vortices outside of the channel are 
not needed to create an effective periodic pinning poten- 
tial. In ratchet channel geometries 1 ^ there are strongly 
pinned vortices outside of the channels; however, we do 
not expect the presence of such vortices to qualitatively 
affect the results we report here. For colloids moving 
through an asymmetric optical trap array, it would be 
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FIG. 6: The vortex velocity (V x ) vs F D for B/B^ = 2.0, 2.25, 
FIG. 5: The average velocity (V x ) vs B/B at F D = 0.03, 2 5j and 2 75 (from top to bottom ). Here the average velocity 
0.025, 0.02, 0.015, 0.01, and 0.005 (from top to bottom). decreases with increasing B. 



possible to remove all colloids outside of the channel so 
that particles are present only inside the ratchet channel. 

To quantify the greater difficulty with which vortices 
flow at higher vortex densities, we examine the average 
normalized vortex velocity (V x ) versus B/B^ for differ- 
ent values of the driving force Fd- In Fig. [T]the initial 
depinning force F c was determined by the initial onset 
of vortex flow; however, the normalized velocity can be 
taken at any value of Fd and in Fig. [5] we plot (V x ) ver- 
sus B/B^, for F D = 0.03, 0.025, 0.02, 0.015, 0.01, and 
0.005. In all cases, (V^) decreases with increasing B/B^ 
expect for some small oscillations caused by commensu- 
ration effects. At Fd — 0.005, (V x ) = for those fields at 
which F c > Fd . The reason that the overall flow velocity 
decreases for increasing B is that the flow patterns orga- 
nize in such a way that only one vortex at a time is able 
to pass through the funnel tip, as will be shown later. 
This is a result of the very large energy cost that would 
be associated with the passage of two or more vortices 
through the tip simultaneously. The single passage con- 
straint causes the flow rate of vortices through the system 
to be roughly constant for fixed Fd , such that for a num- 
ber of moving vortices N m , we find N m cx (V X )B. For 
fixed N m , we obtain (V x ) oc 1/B, which is approximately 
the behavior shown in Fig. [5] The results in Figs. [I] and 
[5] show that jamming phenomena can be realized in su- 
perconducting systems or other systems of particles with 
intermediate to long range interactions. Further, since 
stronger pinning is often a desirable feature for many ap- 
plications of superconductors, some of the concepts from 
studies of jammed systems could be employed to increase 
the effective pinning strength in a device. 



V. DYNAMIC PHASES UP TO B/B^, = 8.0 

We next examine the different vortex dynamical phases 
that can arise in the funnel geometry. In Fig. [5] we plot 




FIG. 7: The vortex positions (dots), funnel geometry (heavy 
lines), and vortex trajectories (light lines) in a small portion 
of the sample at F D = 0.01 for B/B$ = (a) 2.0, (b) 2.25, 
(c) 2.5, (d) 2.75, (e) 3.0, and (f) 3.25. Here the flow at the 
commensurate fields is highly ordered. 



(14) versus F D for B/B^ = 2.0, 2.25, 2.5, and 2.75. The 
depinning force increases with increasing B / B^ over this 
field range, as was seen in Fig. [TJ and above depinning 
(V x ) for a particular value of Fd decreases as more vor- 
tices are added to the system, as was shown in Fig. [5] 

In Fig. [7] we plot the vortex trajectories at fixed 
Fd = 0.01 for B/B^ = 2.0, 2.25, 2.5, 2.75, 3.0, and 3.25. 
The flow is ordered at B/B$ = 2.0 and B/B^ = 3.0 
in Fig. [TJa) and Fig. EJe). The vortices move in fixed 
trajectories and continuously maintain the same neigh- 
bors, indicating that all of the vortices are moving. For 
B/B^ = 2.25, 2.5, and 2.75 in Fig.EKb,c,d), the trajecto- 
ries become increasingly disordered and the vortices no 
longer keep the same neighbors over time, indicative of 
plastic flow. The disordered nature of the flow permits 
the vortices to explore larger regions of phase space, in- 
cluding regimes in which some of the vortices move much 
more slowly than others or become temporarily pinned. 

Fig- Elf) indicates that at B/B^ — 3.25, the vortex 
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FIG. 8: (V x ) vs F D for B/B^ = 3.0, 3.25, 3.5, and 3.75 
(from top right to bottom right), showing the same trends as 
in Fig. 5. 
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FIG. 9: (14) vs F D for (a) B/B^ = 4.0. (b) B/B^ = 4.25. 
(c) B/B^, = 4.5 showing a two step depinning process with a 
cusp near Fd = 0.008 where (14) decreases with increasing 
Fd- (d) B/Bj, — 4.75 showing some small steplike features. 



trajectories are more disordered than at B/B^ = 3.0. In 
Fig. E we plot (14) vs F D for BjB^ = 3.0, 3.25, 3.5, 
and 3.75, which show the same trend observed at and 
above the second matching field. The value of (V x ) at a 
fixed Fd above the depinning threshold decreases with 
increasing B/B^ and the vortex flow becomes partially 
disordered at the incommensurate fields. 

For B I 'B^ > 4.0, the transport becomes more compli- 
cated and transitions between distinct dynamical phases 
begin to occur. The moving phases for B/B^ > 4.0 are 
generally characterized as plastic since a portion of the 
vortices can remain immobile. We find both ordered plas- 
tic motion, where the trajectories of the mobile vortices 
follow a fixed path, and disordered plastic flow phases 
where the vortex motion is more chaotic and the vor- 
tices follow many different paths. In Fig. [Hfa-d) we plot 





FIG. 10: Vortex portions (dots), funnel geometry (heavy 
lines), and vortex trajectories (light lines) in a small portion 
of the sample, (a) B/B,/, = 4.0 at Fd ~ 0.01 where an ordered 
flow occurs, (b) The initial ordered flow at B/ B^ = 4.5 and 
Fd = 0.008, below the cusp in (V x ) in Fig. [9jc), where only 
the vortices along the top wall are moving, (c) The disordered 
flow for BjB^ = 4.5 at Fd = 0.017, above the cusp in (14) in 
Fig. [9{c), where all of the vortices are moving, (d) The high 
drive phase at F D = 0.0275 for B/B^ = 4.5. (e) The initial 
disordered flow at F D = 0.008 for B/B^ = 4.75. (f) The flow 
at F D = 0.01 for B/B^ = 4.75, below the step in (14) which 
occurs near Fd = 0.012 in Fig. EJd). For drives above the 
step in (14) at this field, the vortex flow resembles the flow 
shown in panel (a). 



(14) versus F D for B/B^ = 4.0, 4.25, 4.5, and 4.75. At 
B/Bff, = 4.0, Fig.[9fa) indicates that there is a single well 
defined depinning transition where the vortices flow elas- 
tically, as illustrated in Fig.^UJa). For B/B^ = 4.25, the 
value of F c is depressed as seen in Fig. [9jb) . Here there 
are no sharp jumps in the transport curve, and the vortex 
motion resembles that observed at B/B^ = 4.0 but with 
some additional fluctuations in the vortex trajectories. 

At B/B^ = 4.5, a two step depinning process occurs as 
shown in Fig. [9fc). Initially, only one ordered channel of 
moving vortices forms along one of the funnel walls while 
the remaining vortices are immobile. This is illustrated 
in Fig. fTOT b) for a simulation in which the vortices along 
the upper funnel wall depinned first; depending upon the 
initial random fluctuations, it is also possible for the vor- 
tices along the lower wall to depin first. The remaining 
vortices depin near Fd = 0.008, where a cusp feature ap- 
pears in (V x ). At the cusp, the average vortex velocity 
decreases with increasing Fd, creating a region of neg- 
ative differential conductivity where d(V x )/dFD < 0.0. 
The negative differential conductivity at B/B^ = 4.5 is 
not as pronounced as that observed in simulations and 
experiments with square pinning arraysj^iiii however, we 
find that the type of negative differential conductivity il- 
lustrated in Fig.lHJc) is a common feature at a number of 
the higher order incommensurate fillings we have exam- 
ined in the funnel geometry. In contrast, for the square 
pinning array the negative differential conductivity oc- 
curs only near the first matching field. This indicates 
that it may actually be easier to observe negative differ- 
ential conductivity in a funnel geometry than in a square 
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pinning array. Above the cusp at Fd — 0.008, Fig. [9jc) 
shows that (V x ) smoothly increases with increasing Fd- 

In Fig. [TOTb) we show that the vortex flow below the 
cusp in (V x ) at B/B^ — 4.5 and F D — 0.008 is ordered 
with only a single channel of vortex flow. Figure HOT c) 
indicates that for the same field at Fd = 0.017, above 
the cusp in (V x ), the vortices flow in disordered paths, 
each of which predominantly runs along either the upper 
or lower funnel wall. The transition from only a single 
flowing channel to effectively two flowing channels of vor- 
tices at the second dcpinning transition associated with 
the cusp in (V x ) causes a drop in the mobility of the vor- 
tices due to the competition between the two channels 
for passing a vortex through the bottleneck of the fun- 
nel. Only a single vortex can fit through the bottleneck 
at a time, but vortices in the upper and lower flowing 
channels do not arrive at the bottleneck at synchronized 
times due to the unequal distribution of vortices between 
the two effective channels which causes the two channels 
to flow at different average speeds. As a result, a vortex 
in one channel may reach the bottleneck too soon while 
a vortex from the other channel is still moving through 
the bottleneck, forcing the vortices in the first channel to 
move more slowly until the vortex in the second channel 
has exited the bottleneck and freed it for passage of a 
vortex in the first channel. As Fd increases, the differ- 
ence in flow speed for the two unequally populated chan- 
nels of moving vortices decreases until it is small enough 
that the flow of the two channels becomes synchronized 
on average and a much more orderly passage of vortices 
through the bottleneck occurs, alternating between the 
two channels. In this case, the disordered trajectories be- 
come smoother, as shown for Fo — 0.0275 in Fig. nUT d). 

For B/B^ = 4.75, Fig. [^d) shows that there are sev- 
eral small steps in the (V x ) versus Fd curve which cor- 
respond to changes in the flow. In Fig. ITOT e). the initial 
flow for this field at Fd = 0.008 is partially disordered. 
By Fd = 0.01, shown in Fig- flUT f). a transition to a more 
ordered flow state has occurred where there are two pos- 
sible paths for vortices to follow in the lower channel. At 
higher drives, one of these two lower paths closes, corre- 
sponding to the jump in (V x ) near Fd — 0.012, and the 
flow at higher drives resembles that shown in Fig. ITOT a). 

At B/B,p — 5.0 there is a single depinning transition, 
as illustrated in the plot of (V x ) versus Fn in Fig.[TTJa). 
Right at the depinning transition, the vortices undergo a 
structural transition which is highlighted in Fig. [12] The 
vortex configuration for drives well below depinning con- 
sists of four vortices arranged symmetrically along the 
funnel walls with a fifth vortex in the center of the chan- 
nel near the funnel tip, as shown in Fig. [T2T a). At the 
onset of depinning, illustrated in Fig. [T27b) just below 
Fd = 0.0075, the vortices at the funnel tips move in the 
positive x direction and shift into the adjacent funnel. 
At depinning, which occurs at Fd = 0.0075, Fig. [l~27c) 
indicates the rearrangement that occurs when the center 
vortices move against one of the funnel walls in an alter- 
nating pattern. Above depinning, as shown in Fig. 112( d) 
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FIG. 11: (V x ) vs F D for (a) B/B^ = 5.0. (b) B/B^ = 5.25. 
(c) BjB^ = 5.5. (d) B/B^ = 5.75. 
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FIG. 12: The structural transition that occurs in the single 
stage depinning transition at B/B^, = 5.0 shown in a small 
portion of the sample. Dots: vortex positions; heavy lines: 
funnel geometry; light lines: vortex trajectories, (a) Right 
before depinning, at Fd = 0.006, there are four vortices lining 
the walls and one vortex located in the center of the channel 
near the funnel tip. (b) The first rearrangement just below 
Fd = 0.0075 occurs when the center vortices move through 
the aperture in the positive x direction and pass into the 
adjacent funnels, (c) At depinning, which occurs at Fd = 
0.0075, these center vortices move against the funnel walls in 
an alternating pattern associated with a transient flow, (d) 
Above depinning at Fd = 0.012, the vortices flow strictly 
along the funnel walls. 



for Fd = 0.012, the vortices flow strictly along the funnel 
walls in two channels that never mix. 

At B/Btp = 5.25, 5.5, and 5.75, shown in Fig.rjjTh), 
fTTTcj . andfTTTd). respectively, there can be multiple dy- 
namical transitions between ordered and disordered flow 
phases. At B/B^ = 5.25, the flow at large Fd is very 
similar to that shown for B/B^ = 5.0 in Fig.[T2Td): how- 
ever, the flow initiates at a much lower drive and takes 
the form of an ordered winding channel in which about 
40% of the vortices are moving, as shown in Fig. ITST a). 
As the vortices in the winding channel move through the 
system, the vortices pinned along the funnel walls un- 
dergo an oscillatory motion which is most easily seen in 
Fig. null a) for the pinned vortex near the center of each 
long funnel wall. This oscillation is a response to the pas- 
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FIG. 13: The vortex positions (dots), funnel geometry (heavy 
lines), and vortex trajectories (light lines) in a small portion 
of the sample, (a) The initial flow at Fd = 0.006 for B/B$ = 
5.25. (b) The initial depinning at F D = 0.012 for B/B^ = 
5.75. (c) The trajectories in the fluctuating phase for Fd = 
0.019 at B/Btf, = 5.75. (d) The high drive ordered phase at 
F D = 0.032 for B/B^ = 5.75. 
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FIG. 14: (V x ) vs F D for (a) B/B^ = 6.0. (b) B/B^ = 6.75. 



sage of an individual vortex through the flowing channel; 
the pinned vortex shifts in the positive x direction as the 
moving vortex approaches, and shifts back again after the 
moving vortex has passed. The vortices at the corners of 
the funnel undergo little to no shift. In Fig. UTT d) at 
B/B^ = 5.75 the transitions between different dynami- 
cal phases are associated with changes in the fluctuations 
of the velocity signal. At this filling, the sharp depinning 
transition takes the system into the ordered braiding flow 
phase illustrated in Fig. HSJb), a state with low levels of 
velocity fluctuations. A transition to a highly fluctuating 
phase occurs near Fd = 0.15 and corresponds to the on- 
set of the disordered flow phase illustrated in Fig. [T"3tc). 
For higher drives, a transition to an ordered phase occurs 
near Fd = 0.0225 when the trajectories become partially 
ordered. A velocity jump near Fd = 0.025 marks the 
transition to the completely ordered flow phase shown in 
Fig.HSd). 

At B I B,/, — 6.0, where a commensuration peak in F c 
was missing in Fig. [TJ the velocity-force curve shows a 
single depinning transition into a random flow phase, as 
shown in Fig. [Uf a) . No sharp transitions between differ- 
ent phases appear; however, the general form of the vor- 
tex flow changes gradually as Fd increases. Figure [TBT a) 
shows an initial shift in vortex positions occurring at a 
drive just below the depinning transition. The vortex 
at the tip of the leftmost funnel has pushed the vortex 
at the base of the adjacent funnel into a position along 
the lower wall in a pattern that repeats every two funnel 




FIG. 15: The vortex positions (dots), funnel geometry (heavy 
lines), and vortex trajectories (light lines) in a small portion 
of the sample, (a) The initial motion at Fd = 0.006 for 
B/B,f, = 6.0. (b) The flow above depinning at F D = 0.008 
for B/B^ = 6.0. (c) The flow at F D = 0.019 for B/B^ = 6.0. 
(d) The flow at Fd = 0.028 for B/Bj, = 6.0. (e) The ordered 
vortex flow at Fd = 0.025 for BjB^ = 6.75, below the cusp 
in (V x ) vs Fd seen in Fig. |14f b). (f) The disordered flow at 
F D = 0.028 for BjB^, = 6.75, just after the cusp in (V x ). 



plaquettes, indicating that the vortices in the center of 
the channel were not well pinned, resulting in lack of a 
peak in F c at this field. Just above depinning, shown 
at Fd = 0.008 in Fig. HBT b). a disordered flow occurs in 
which the two vortices at the corners of each funnel do 
not participate. Although the flow is disordered, there 
are clearly defined regions inside the funnels which the 
vortices completely avoid, as shown by the lack of tra- 
jectories passing through large areas of the funnels. As 
Fd is further increased, the trajectories become increas- 
ingly disordered and the vortices in the corners of the 
funnels begin to take part in the motion, as shown in 

Figs.niM). 

More clearly defined transitions between ordered and 
disordered flow states occur at B/B^ = 6.25, 6.5, and 
6.75, as illustrated in Fig. \W[b) for B/B,p = 6.75. 
The transitions are characterized by cusp structures in 
(V x ) associated with negative differential conductivity. 
Fig- UST e) shows the initial ordered flow at Fd = 0.028 
for B I B^ = 6.75, below a pronounced cusp in (V x ). Here 
the two vortices in the corners of each funnel remain im- 
mobile so that the flow is actually plastic. These vortices 
become mobile at the transition to the disordered phase, 
illustrated at Fd — 0.028 in Fig. [TSTe). A similar set of 
dynamics appears at B/B^ = 6.25 and B/B^ — 6.5. For 
driving forces higher than those we have examined, it is 
possible that further dynamical transitions could occur 
at which different ordered phases arise. 

Depinning at B/B^ = 7.0 occurs in a single step, 
as illustrated in Fig. [TBT a) where we plot (V x ) versus 
Fd- At this field, Fig. [JJ indicates that there is a peak 
in F c . Above depinning, the vortices flow in the or- 
dered pattern shown in Fig. 117( a). where two vortices 
remain pinned in the corners of each funnel. For drives 
Fd > 0.3 higher than those we consider here, a transi- 
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FIG. 16: (V x ) vs F D for (a) B/Bj, = 7.0 and (b) B/B4, 
7.75, where three distinct moving phases appear. 
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FIG. 17: The vortex positions (dots), funnel geometry (heavy 
lines), and vortex trajectories (light lines) in a small portion 
of the sample, (a) The vortex flow state at Fd = 0.014 for 
B/B4, = 7.0. (b) The braided flow phase at F D = 0.008 for 
BjBcf, = 7.75. (c) The second flow phase at Fd = 0.011 for 
B/B4, = 7.75. (d) The third flow phase at F D = 0.021 for 
B/B^ = 7.75. 



tion to a disordered flow state may occur once the drive 
is large enough to cause the two immobile vortices to de- 
pin. At B/Bcf, = 7.25, 7.5, and 7.75, multiple steps occur 
in the velocity-force curves, as illustrated in Fig. fToTb) 
for B/B^ = 7.75. Each step is associated with a dis- 
tinct type of ordered flow. At low drives, we find the 
braided flow shown in Fig. HTJb) at F D = 0.008. Above 
the first low step in (V x ) near Fd = 0.01, an alternating 
braided flow occurs in which every other plaquette con- 
tains two possible flow paths along the lower funnel wall, 
as illustrated in Fig.fTTtc) at Fo = 0.011. Above the sec- 
ond larger step in (V x ) near Fd = 0.012, the completely 
ordered flow phase shown in Fig. fTTTd) at Fd = 0.021 
occurs. 



VI. DYNAMICS AT HIGHER FIELDS B/B^ > 8 

For higher fields B/B^ > 8.0 we observe the same 
trends found for fields just below B/B^ = 8.0. The in- 
teger matching fields typically display a single ordered 
flow phase, while at the incommensurate fields, multi- 
ple flow phases occur with ordered-ordered or ordered- 
disordered flow transitions. Another trend is that as the 
field increases, the number of immobile vortices appear- 
ing in the ordered flow phases increases. In Fig. [T5] we 
plot (V x ) versus F D for B/B^ = 8.0, 8.25, 8.5, and 8.75, 
where the trend outlined above can be seen. There is a 
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FIG. 18: (V x ) vs F D for (a) B/B^ = 8.0, (b) B/B^ = 8.25, 
(c) B/B^ = 8.5, and (d) B/B^ = 8.75. 
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FIG. 19: The vortex positions (dots), funnel geometry (heavy 
lines), and vortex trajectories (light lines) in a small portion 
of the sample, (a) The ordered flow phase at Fd = 0.014 
for B/B$ = 8.0. (b) The disordered flow at F D = 0.017 for 
B/B,f, = 8.5. (c) The initial partially ordered flow phase at 
F D = 0.008 for B/B^ = 8.75. (d) The disordered flow phase 
with two immobile vortices in the corners of each funnel at 
F D = 0.025 for B/B^ = 8.75. 



single depinning step in Fig. [TSJa) at B/B^ = 8.0 into 
the ordered flow phase shown in Fig.fTWa) at Fd = 0.014 
where two vortices remain pinned at the corners of each 
funnel. Figure [TW b) shows the disordered flow at the 
higher drive F D = 0.017 for B/B^ = 8.5. At this field, 
there is a transition near Fd — 0.025 to an ordered phase 
with vortex trajectories that are very similar to those 
shown in Fig. ITWa) for B/B,/, — 8.0. The initial flow at 
B I Bfj, — 8.75 is partially ordered with four immobile vor- 
tices in each funnel plaquette, as shown at Fd — 0.008 
in Fig. [TWc). At higher drives for this field, the flow 
becomes disordered but only two of the immobile vor- 
tices begin to move, leaving two vortices immobile in 
the corners of each funnel, as shown in Fig. HW d) for 
F D = 0.025. 

We find a distinctive integer matching field at which 
there is not merely a single dynamic flow phase state. 
This field, B/B^ — 9.0, also fails to produce a depinning 
peak as shown in Fig. [TJ We plot (V x ) versus Fd for 
B/B^ = 9.0 in Fig. I2"07 a). which indicates that there are 
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FIG. 20: (14) vs Fd for (a) B/B^ = 9.0 and (b) B/B^ = 
9.75. 




FIG. 21: The vortex positions (dots), funnel geometry (heavy 
lines), and vortex trajectories (light lines) in a small portion 
of the sample at B/Bj, = 9.0 showing braided flows, (a) First 
flow regime shown at Fd = 0.01. (b) Second flow regime 
shown at Fd = 0.014. (c) Third flow regime shown at Fd = 
0.023. 



three distinct moving phases. All three phases involve 
braided vortex flows, but the details of the braiding and 
the number of vortices participating in the flow varies 
as a function of drive, as shown in Fig. [2T] where the 
three different braided flows are illustrated. The sup- 
pression of the depinning peak at B/B^ = 9.0 is related 
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FIG. 22: The vortex positions (dots), funnel geometry (heavy 
lines), and vortex trajectories (light lines) in a small portion 
of the sample, (a) The initial motion at Fd = 0.008 for 
B/Bj, = 9.75. (b) The second flow phase at F D = 0.012 for 
B/B,p = 9.75. (c) The higher drive ordered flow phase at 
F D = 0.028 for = 9.75. (d) The flow at F D = 0.017 for 

BjB$ = 11.0. (e) A complex partially ordered flow pattern 
at Fd = 0.017 for B/B^ = 11.25. (f) A higher drive ordered 
alternating flow pattern at Fd = 0.025 for B/B^ = 11.75. 



to the suppression of the peak at B/B^ = 6.0. At each 
of these fields, the system exhibits an instability in which 
the number of vortices present in the center of the chan- 
nel for zero applied drive is different than the number 
of vortices that can be stabilized along the center of the 
channel just above depinning. At B/B^ = 6, in each 
funnel plaquette there are two vortices positioned along 
the center line of the channel at zero drive, but under fi- 
nite drive one of these vortices moves against the funnel 
wall, leaving only one vortex in the center of the channel. 
This transition can be seen by comparing Fig. 2(e) with 
Fig. 15(a,b). At this filling, there is only a small energy 
difference between placing two vortices at the center of 
the channel and placing only one vortex at the center 
of the channel with the other vortex against the funnel 
wall. This produces an instability which contributes to 
the disordered flow we find at B/B^ — 6. An opposite 
transition occurs for B/B^ = 9. As seen in Fig. 3(b) 
and Fig. 21(a), there is only one vortex per plaquette 
in the center of the channel at zero applied drive, but 
above depinning a second vortex moves into the center 
of the channel and is stabilized there. The existence of 
this change in the occupancy of the center of the channel 
at depinning from one to two or two to one is associ- 
ated with the suppression of commcnsuration peaks in 
the critical depinning force at the fields B/B^ — 6 and 
BjB^ = 9 seen in Fig. 1. 

A series of flow transitions also occurs at B/B^ = 9.75 
in Fig. 1217T b). corresponding to the dynamical flows il- 
lustrated in Fig. l22T a.b.c). In the initial flow phase at 
low drives, Fig. [22T a) indicates that only a portion of 
the vortices are moving in a path that passes through 
the bottom half of each funnel, forming a single flowing 
channel. At higher drives a transition to the symmetric 
flow state illustrated in Fig. 122T b) occurs. At still higher 
drives, an intermediate disordered phase appears before 
the flow reorders into the state shown for Fjj = 0.028 in 
Fig. I22Tc). where there are two immobile vortices on the 
bottom half of each funnel and one immobile vortex on 
the top half of the funnel, with flow occurring through 
two well-defined channels. 

At higher fields, similar dynamical regimes occur and 
larger numbers of vortices can become immobile along 
the funnel walls, such as at B/B<p = 10 where there 
is a single step depinning transition into a state with 
four immobile vortices in each funnel. It is also possible 
for vortices to be immobilized in the center of the chan- 
nel rather than along the funnel walls, as illustrated for 
B/B^ = 11.0 in Fig. !22Td). At this field, a single step 
depinning transition occurs into the state shown where a 
total of three vortices are immobilized: two in the corners 
of the funnel and one right in the funnel center. As B/B^ 
increases, additional partially ordered phases occur with 
increasingly intricate flow patterns such as the one shown 
in Fig.[22Te) at B/B^ = 11.25. We also find phases with 
asymmetric flow where the asymmetry alternates from 
one funnel to the next, such as the state illustrated in 
Fig. |22T f) for B/B^ = 11.75. Here the leftmost funnel 



13 




'le-05 



0.0001 



0.001 



0.01 



FIG. 23: The power spectrum S(v) versus frequency v in 
inverse simulation time steps. Lower curve: The ordered 
phase at Fd = 0.02 for B/B^, = 5.0, where a series of 
peaks occurs due to the periodic velocity signal of the vor- 
tices. Upper curve: The disordered flow phase at Fd = 0.02 
for BjB^, = 6.0, where a broadband noise signal arises with 
a 1/u 2 ' 5 form. The same types of spectra appear for ordered 
and disordered flow phases at other fields. 



contains two flowing channels of vortices interacting with 
the upper funnel wall and one channel interacting with 
the lower funnel wall, while in the middle funnel this flow 
pattern is reversed. 

We expect the trends described above to continue for 
higher matching fields. Eventually, however, there may 
be a crossover to a state where the vortex lattice con- 
stant is smaller than the width of the funnel aperture, 
at which point it would be possible to move two vor- 
tices through the funnel tip at the same time without 
much of an additional energy cost. In this case, the av- 
erage velocity might sharply increase since the flow in 
the aperture would transition from one-dimensional to 
quasi-two-dimensional, and it may be possible that fur- 
ther oscillations in the average velocity would arise at 
even higher fields as larger groups of vortices could pass 
through the aperture simultaneously for increasing mag- 
netic field values. 



VII. NOISE MEASUREMENT AND 
TEMPERATURE EFFECTS 

With current experimental noise measurement tech- 
niques, it should be possible to characterize the differ- 
ence between ordered and disordered flow phases. The 
vortex velocities we measure in our simulations are pro- 
portional to the experimentally measured vortex voltage 
signal, which can be analyzed using the power spectrum 
defined as 



S(u) = 



1 

2tt 



V(t)e 



'dt 



(3) 



In Fig.[23]we plot S(v) for B/B^ = 5.0 at F D = 0.02 in a 
regime of ordered flow. We find narrow band noise peaks 
at the characteristic frequencies of the ordered oscillatory 
motion of the vortices. Also plotted in Fig. [231 is S(u) for 
BjB^ = 6.0 at Fd = 0.02, where the vortex flow is dis- 
ordered. Here the power spectrum has a broad band or 
1// Q noise characteristic with a = 2.5. In general, for 
other values of B/B$, the ordered phases produce nar- 
row band spectra while the disordered phases produce 
broad band noise signals. The appearance of the narrow 
band noise in the ordered phase implies that phase lock- 
ing phenomena could be induced by adding an ac drive 
component to the dc driving force. When the frequency 
of the ac drive matches the intrinsic or higher harmonic 
frequencies of the moving vortices in the ordered phase, a 
series of steps should appear in the velocity force curves. 

We next consider the effect of temperature. In pe- 
riodic pinning arrays, the most pronounced matching 
field peaks observed in experiments occur at tempera- 
tures near T c . This has been attributed to the thermal 
suppression of the intrinsic pinning in the sample that 
competes with the periodic pinning. Another factor that 
could play a role is that the vacancies or interstitials 
present in the vortex lattice away from commensurate 
fields can be relatively mobile under thermal activation, 
and so increasing the temperature could depress the de- 
pinning force to a larger degree at incommensurate fields 
than at commensurate fields where vacancies or intersti- 
tials are not present. We find that for a temperature of 
T = 0.1, which is well below the melting temperature of 
the vortex lattice, the depinning force peaks at the com- 
mensurate fields are robust while the depinning force at 
the incommensurate fields is depressed. The peaks in F c 
that appear at B/B^ = 4.5 and B/B^ = 5.5 for T = 
are absent at T = 0.1, indicating that for finite tem- 
peratures, pronounced peaks are only observable at the 
matching fields whereas submatching field peaks are less 
robust. The general trend of increasing F c with increas- 
ing B/Bif, also remains robust at finite temperature. 



VIII. DISCUSSION 

Although our work is focused on a superconducting 
vortex system, we expect our results to be general for 
other systems of particles with repulsive interactions. For 
example, in colloidal systems where the volume density 
is sufficiently low that steric contact of the colloidal par- 
ticles does not occur, many of the same results should 
still apply. One aspect of the vortex system that we do 
not take into account is the possibility of multi-quanta 
vortex formation under certain conditions. In static sys- 
tems, the formation of giant vortices has been observed in 
cases where the giant vortices produce a higher symmetry 
of the vortex configuration, such as in quasiperiodic pin- 
ning arrays^ It is possible that dynamical multi-quanta 
vortices could form in a funnel geometry. For example, 
if two vortices are forced into close proximity near the 
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funnel tip, they could merge to form a two-quanta vor- 
tex which would then pass through the funnel tip. The 
formation of such dynamical multi-quanta vortices could 
produce interesting signatures in the transport curves. 
The experimental system closest to the system studied 
here is periodic asymmetric channel geometries- - where 
small oscillations in the critical current due to commen- 
suration effects were observed. Since these commensura- 
tion effects were very weak, it is possible that they were 
due to edge effects. More recent experiments with ge- 
ometries that avoid the effect of edges have now revealed 
much more prominent commensuration effects ^ indicat- 
ing that some of the phases we observe may be occurring 
in this system. 

IX. SUMMARY 

In summary, we have used numerical simulations to 
examine the vortex configurations and dynamics in a pe- 
riodic funnel array. The vortex configurations we observe 
are generally different from those found for a single iso- 
lated triangular sample due to the coupling between vor- 
tices in adjacent funnels. As a function of field we find a 
series of depinning threshold peaks at the matching fields 
where the vortex configurations are ordered. In some 
cases, matching peaks are missing due to the fact that 
the vortex configuration contains pairs of vortices that 
are located close together near the funnel apertures. We 
also observe a general increase in the depinning thresh- 
old with increasing vortex density, which is opposite from 
the normal trend for vortices in two-dimensional pinning 
arrays where the vortices are directly trapped by pin- 
ning sites. In the funnel geometry, the pinning is a result 
of the repulsive vortex- vortex interaction forces, and as 



the vortex density increases it becomes more difficult for 
the vortices to overcome these repulsive forces and flow 
through the funnel tip. We also find a rich variety of 
dynamical phases, including ordered elastic and ordered 
plastic phases where the vortices follow fixed trajecto- 
ries and disordered phases where vortices mix chaotically. 
The phases generally organize in such a way that only 
one vortex passes through the tip of the funnel at a time. 
Due to this constraint, the average velocity of an indi- 
vidual vortex decreases with increasing field such that 
the sum of the velocities of all of the vortices at fixed 
drive remains close to constant for increasing field rather 
than increasing with increasing field. This behavior is 
similar to the response of grains in an hourglass. Tran- 
sitions between the different dynamical phases appear as 
jumps or cusps in the velocity-force curves, and there are 
even regimes where the average vortex motion decreases 
with increasing drive, producing a negative differential 
conductivity. At higher fields, moving states can form in 
which a single vortex remains immobile at the center of 
a funnel while other vortices flow around it. In general, 
ordered flow phases occur at the matching fields, while 
at non-matching fields the flow is disordered for at least 
some regime of driving forces. The ordered phases are 
associated with sharp narrow band or washboard veloc- 
ity noise signals while the disordered phases have l/f a 
velocity noise spectra. Our results should also be gener- 
alizable to other systems of repulsively interacting parti- 
cles moving through a funnel geometry, such as colloids 
or Wigner crystals. 

We thank B. Plourde for helpful discussions. This 
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the U.S. DoE at LANL under Contract No. DE-AC52- 
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